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ABSTRACT 

The Palomar Integral Field Spectrograph was used to probe a variety of envi- 
ronments in nine nearby galaxies that span a range of morphological types, lumi- 
nosities, metallicities, and infrared-to-blue ratios. For the first time, near-infrared 
spectroscopy was obtained for nuclear or bright H II regions in star-forming 
galaxies over two spatial dimensions (5'/7 x lO'.'O) in the [Fe II] (1.257 jttm), 
[Fe II] (1.644 /jm), Pa/3 (1.282 /mi), H 2 (2.122 /mi), and Br7 (2.166 /mi) transition 
lines. These data yield constraints on various characteristics of the star-forming 
episodes in these regions, including their strength, maturity, spatial variability, 
and extinction. The H II regions stand out from the nuclei. Unlike observations 
of nuclear regions, H II region near-infrared observations do not show a spatial 
coincidence of the line and continuum emission; the continuum and line maps of 
H II regions usually show distinct and sometimes spatially-separated morpholo- 
gies. Gauging from Pa/3 and Br7 equivalent widths and luminosities, the H II 
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regions have younger episodes of star formation than the nuclei and more intense 
radiation fields. Near-infrared line ratio diagnostics suggest that H II regions 
have "purer" starbursting properties. The correlation between ionizing photon 
density and mid-infrared color is consistent with the star formation activity level 
being higher for H II regions than for nuclei. And though the interpretation 
is complicated, on a purely empirical basis the H II regions show lower Fe 1+ 
abundances than nuclei by an order of magnitude. 

Subject headings: galaxies: ISM — infrared: galaxies — galaxies: individual 
(IC 10, NGC 693, UGC 2855, NGC 1569, NGC 2388, NGC 4418, NGC 6946, 
NGC 7218, NGC 7771) 

1. INTRODUCTION 

Stellar processing is the main evolutionary process affecting the cosmos. Nucleosynthesis 
and mass return to the interstellar medium modify the chemical composition of the Universe 
and the fraction of mass bound in stars. Star formation dominates the present-day radiation 
field, and most of it happens in normal galaxies (e.g. Kim & Sanders 1998; Driver 1999). 
Relatively little, however, is known about star formation on the scale of a galaxy, including 
its drivers and inhibitors, and the role of interactions (e.g. Volker 2000; Dopita et al. 2002; 
Hameed & Young 2003). The ISO Key Project on the Interstellar Medium of Normal 
Galaxies (Helou et al. 1996; Dale et al. 2000) aimed to better understand the large-scale 
physics of the interstellar medium through an array of mid- and far-infrared data on a diverse 
sample of 69 nearby star-forming galaxies (e.g. Dale et al. 1999, Dale et al. 2000, Helou et 
al. 2000, Hunter et al. 2001, Malhotra et al. 2001, Contursi et al. 2002, Lu et al. 2003). In 
this contribution we present and discuss near-infrared integral field spectroscopy for several 
of these objects. 

For optical depths A v > 3, the brightest emission lines that probe the star- forming 
interstellar medium are found in the near-infrared. Vibrationally excited H 2 emission at 
2.122 fim arises by collisional excitation from shocks or radiative excitation in intense ultra- 
violet radiation fields (e.g. Goldader et al. 1997), or possibly from reprocessed X-ray illumi- 
nation (Maloney, Hollenbach, & Tielens 1996). On the other hand, the [Fe II] lines at 1.257 
or 1.644 /im probe supernova-shocked gas (the final stage of processing) and hard power- 
law/X-ray environments, both of which can release iron atoms and ions through interstellar 
dust grain sputtering and evaporation; recent HST work on M 82 and NGC 253 by Alonso- 
Herrero et al. (2003) indicates that as much as 70% of the [Fe Il]1.644 /im flux ultimately 
derives from supernovae, and only 6-8% from H II regions. In contrast to [Fe II], hydrogen 
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lines like Pa/3 (1.282 /jm) and Br7 (2.166 fim) directly trace the gas ionized by young mas- 
sive stars. Thus, an [Fe II]-to-hydrogen line ratio can be interpreted as an indicator of the 
maturity of the local star formation episode, as reflected in the ratio of supernova-shocked 
gas to molecular gas feeding star formation. Moreover, the ratios [Fe II] 1.257 /xm/Pa/3 and 
H2(2.122 /im)/Br7 discriminate between shock excitation from supernova remnants or hard 
X-ray heating from power-law sources (Larkin et al. 1998), and are essentially unaffected by 
reddening. 

In addition to studying these diagnostics, we use the near-infrared line and continuum 
images to explore the maturity and spatial progression in the star formation. Coupled 
with Ha imaging, the near-infrared hydrogen line fluxes are used to estimate the optical 
depth and reconstruct the intrinsic emission, total ionizing flux, and other properties of the 
local interstellar medium. Data at other wavelengths help to further probe the interstellar 
medium. For example, new optical spectroscopy and archival mid-infrared imaging allow us 
to investigate trends with metallicity and infrared colors. The various results will point to 
important physical differences between nuclear and extranuclear/H II regions in star-forming 
galaxies. 

2. THE SAMPLE 

The targets (Tables 1 and 2 and Figures 1-5) derive from the ISO Key Project on 
the Interstellar Medium of Normal Galaxies. The Key Project collected infrared data for 
69 star-forming galaxies: carbon, oxygen, and nitrogen fine-structure line fluxes between 50 
and 160 /xm (Malhotra et al. 2001), mid-infrared spectra between 3 and 12 /im (Lu et al. 
2003), and mid-infrared maps at 7 and 15 /im (Dale et al. 2000). In addition, the following 
ancillary ground-based data have been obtained for a large portion of the sample: broadband 
B and narrowband Ha imaging, long-slit optical spectroscopy, and literature CO fluxes. 

The Key Project sample had the original following criteria: /^(60 /im) > 3 Jy; a 
published redshift; and no AGN or Seyfert classification on NED (at that time). The sample 
explores the full range of morphology, luminosity (Lfir from less than 10 8 L Q to as large 
as 1O 12 L ), infrared-to-blue ratio (0.05 to 50) and infrared colors (see § 2 of Dale et al. 
2001 for a complete description of the sample selection). The subset of the Key Project 
sample selected for this particular ground-based follow-up project was chosen according to 
three criteria. The subset should: contain targets of high 7 /im surface brightness (to ensure 
detectable targets in the near- infrared); span a range of properties like metallicity, infrared- 
to-blue ratio, mid-infrared color, morphology, etc.; and be constrained in redshift such that 
multiple galaxies could be observed with the same grating tilts, to minimize observational 
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overheads. 



2.1. 



IC 10 



IC 10 is classified as a Magellanic-type dwarf irregular galaxy and has been extensively 
studied in the optical (e.g. Lequeux et al. 1994; Hunter 2001). The galaxy appears to exhibit 
ongoing star formation at a rate of 0.15 M yr -1 (Thronson et al. 1990) and has the highest 
global surface density of Wolf-Rayet stars in the Local Group (Massey & Armandroff 1995). 
The H I emission spans some seven times the optical extent of the galaxy and Wilcots 
& Miller (1998) suggest that IC 10 is still forming through accretion of the surrounding 
material. 



NGC 693 is a non-Magellanic irregular system with average infrared-to-optical and 
infrared colors. This galaxy's fairly high global / 1/ (6.75/xm)// 1/ (15/xm) ratio is presumably 
indicative of a system dominated by diffuse cirrus emission (a low / J/ (6.75/xm)// I/ (15/xm) ratio 
would indicate intense star formation; Dale et al. 2000). This impression is supported by the 
optical spectrum (integrated along the major axis), which shows strong stellar absorption 
features in the blue (e.g. calcium H and K, and higher-order Balmer lines). 



UGC 2855 is a late-type spiral with prominent arms that are quite striking in the mid- 
infrared (Dale et al. 2000). Hiittemeister, Aalto, & Wall (1999) argue that this galaxy lacks 
evidence for strong shocks and is in a pre-burst stage. Their claim is based on the ~ 8 kpc- 
long continuous molecular bar, which is relatively undisturbed and rotating in a solid-body 
fashion. Moreover, the 12 CO(l — > 0)/ 13 CO(1 — > 0) ratio is low and the Ha emission along 
the bar is weak. Our optical spectrum of a representative H II region in UGC 2855 shows 
very strong nebular emission lines superposed on a weak and relatively featureless underlying 
continuum. UGC 2866 is a companion galaxy. 



2.2. 



NGC 0693 



2.3. 



UGC 2855 
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2.4. NGC 1569 

The post-starburst NGC 1569 is famous for its two brightest compact super-star clusters 
("A" and "B"). In this study we have targeted two of the galaxy's H II regions. The 
southeastern and northwestern H II regions respectively correspond to regions 7 and 2 of 
Waller (1991); the latter has also been designated as "Region C" by Greve et al. (1996) 
and "GMC 3" by Taylor et al. (1999). NGC 1569 has unusually strong ionic line intensities 
relative to measures of the dust and photodissociation region emission, in both the mid- and 
far-infrared. It stands out as extreme among the general Key Project sample in the ratios 
of [O III] 88 //m/far-infrared, [O III] 88 //m/[C II] 158 /im, and [O III] 88 //m/[0 I] 63 //m, 
consistent with the idea that the galaxy is powered by vigorous ongoing star formation 
(Malhotra et al. 2001). The extreme mid-infrared colors point to an intense radiation field 
(Dale et al. 2000). The 2.5 — 11.6 /im ISOPHOT spectra from 24"x24" apertures centered on 
the southeastern and northwestern H II regions (but overlapping with the super-star clusters) 
are dominated by extraordinarily strong ionic emission, in this case the [S IV] 10.5 /im line, 
relative to the underlying dust continuum emission. The optical spectra of H II regions in 
this galaxy display extremely strong [O III] 4959, 5007 A emission due to the high electron 
temperature T e , a result of the low abundance of oxygen, the primary coolant. However, since 
the infrared lines (see above) are not sensitive to T e , the strong line-to-continuum contrast 
contrast in the infrared must have a different origin, possibly a deficiency of interstellar dust 
emitting at these wavelengths due to the extremely intense ultraviolet radiation field (see 
Lu et al. 2003). 



2.5. NGC 2388 

NGC 2388 is a spiral galaxy with fairly unremarkable infrared properties. We did not 
obtain an optical spectrum of NGC 2388. 



2.6. NGC 4418 

NGC 4418 has the most unusual infrared properties in our sample. Roche et al. (1986) 
interpret the deep 9.7 /im spectral trough in its spectrum as a silicate absorption feature, and 
infer a high extinction, Ay > 50, towards the nuclear source. The galaxy has an extremely 
small / I/ (6.75/im)// I/ (15 / um) ratio and an extremely large J^ 6 ^™^ ratio, indicative of an 
intense interstellar radiation field (Dale et al. 2000). Its mid-infrared spectrum lacks the 
typical polycyclic aromatic hydrocarbon features (Lu et al. 2003). NGC 4418 also exhibits 
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an unusually low [C II] /far- infrared ratio (an upper limit, in fact); the interpretation favored 
by Malhotra et al. (2001) is that interstellar dust grains are highly positively charged in 
this active galaxy, leading to less efficient heating by photo-ejected electrons from dust 
grains. The optical spectrum of NGC 4418 is dominated by spectral features from late-type 
stars. The only emission line evident is a weak emission peak near 6620 A which may be 
a component of a complete Ha profile that includes a small absorption dip at the galaxy's 
systemic velocity. Due to the large extinction towards the infrared-active nuclear region of 
this galaxy, the optical spectrum does not probe this region, but instead simply presents a 
spectrum characteristic of the older stellar population of the galaxy's disk. 

2.7. NGC 6946 

NGC 6946 is a well-studied, nearly face-on galaxy with classic spiral structure and a 
starbursting core (Engelbracht et al. 1996). Schlegel, Blair, & Fesen (2000) find 14 pointlike 
X-ray sources in NGC 6946, the brightest of which is unremarkable in the mid-infrared 
(compare with the mid-infrared maps in Dale et al. 1999). The H II region we targeted 
to the NE shows up in the ROSAT HRI image as a faint arc ~0f5 in length, whereas the 
nucleus is associated with a point source with the third highest X-ray luminosity in the 
galaxy (L x (0.5 - 2.0 keV) ~ 29 • 10 30 J). 

2.8. NGC 7218 

NGC 7218 is a late-type spiral with typical infrared properties. The galaxy has a short 
bar oriented roughly perpendicular to the major axis. The optical emission lines show a level 
of star-forming activity intermediate between NGC 693 and NGC 7771 on the one hand, and 
UGC 2855 and NGC 1569 on the other. 



2.9. NGC 7771 

NGC 7771 has a prominent circumnuclear star-forming ring of major axis diameter 6- 
7" with psIO radio-bright clumps and ~10 near-infrared-bright clumps (Smith et al. 1999; 
Reunanen et al. 2000). The radio and near-infrared sources in the ring are offset from each 
other, lending to the idea that there are multiple star-forming populations of different ages. 
A possible interaction with the nearby dwarf NGC 7770 may be related to much of this 
activity. Smith et al. (1999) find that the radio regions are ~4 Myr old, whereas the near- 
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infrared spots are slightly older, some psIO Myr. Reunanen et al. (2000) present Br7 and 
H 2 line maps and find that the starburst ring exhibits small Br7 equivalent widths (of order 
unity dex) and their interpretation is an instantaneous starburst that occurred 6-7 Myr ago. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. Near-Infrared Spectroscopy 

The spectroscopic observations were carried out between September 2000 and Novem- 
ber 2001 (Table 2) using the Palomar Integral Field Spectrograph (PIFS; Murphy, Matthews, 
& Soifer 1999) on the Palomar 200-inch Telescope. PIFS provides a 5'/4x9 / . / 6 field of view 
using eight separate / . / 67x9'/6 slits to feed two independent spectrographs within the same 
liquid N 2 -cooled dewar. Each slit is four pixels wide and 58 pixels long (0'.'167 pixel -1 ). The 
high resolution mode (R ~ 1000 — 1500) was used, providing a usable wavelength range 
of 0.03-0.04 /im and a velocity resolution of Av ~ 200 — 300 km s -1 near the observed 
wavelengths. 

The standard observational procedure consisted of a set of four five-minute on-off inte- 
grations. Thus 20 minutes were spent integrating on source, and an equal amount of time 
off source. The only exceptions to this procedure involved UGC 2855 H 2 and NGC 693 Br7, 
for which we observed 40 minutes and 45 minutes on source, respectively. For all observa- 
tions, the b"A^'.Q field of view was oriented with the long axis at a position angle of 90°. A 
positional dithering pattern was employed for the sequence of observations to recover seeing- 
limited spatial sampling in the cross-slit direction. This dithering is the reason the resultant 
maps (5'/7x lO'.'O) are slightly larger than the PIFS field of view. Spectral calibration lamp 
data were taken immediately after each set of four dithered integrations. Corrections for 
atmospheric opacity and spectral flat-fielding are derived from observations of F8-G9 main 
sequence stars. An estimate of the near-infrared PSF was periodically obtained, with each 
measurement comprised of four dithered exposures of a nearby field star using the PIFS 
imaging mode (see Table 2). 

The atmospheric conditions for each set of observations are listed in Table 2. Since about 
half of the observations were taken in non-photometric conditions, for uniformity all the 
observations are flux-calibrated using 2MASS JHK data, via fluxes from cutouts matched 
to the PIFS field of view. This technique is uncertain. First, near-infrared broadband fluxes 
are contaminated at some level by nebular line emission. In very young starbursts/blue 
compact galaxies, for example, the line contribution can be as large as ~5-10% of the total 
near-infrared flux (e.g. Figure 3 of Kriiger, Fritze-v. Alvensleben, & Loose 1995; Vanzi et 
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al. 2000). However, except for IC 10 Pap, NGC 1569 NW Pa/3, and NGC 6946 NE Br 7 , 
none of the equivalent widths for our sample indicate broadband flux contributions larger 
than 2%. Second, the 2MASS and PIFS aperture positionings may not exactly coincide. 
Third, depending on the location of the emission within the PIFS field of view, 2MASS 
flux dilution can be a factor since the 2MASS images are ~ 3" in resolution whereas the 
PIFS data typically are sub-arcsecond. Based on a comparison of 2MASS flux calibration 
and the calibration we obtained from standard stars on the photometric nights, we estimate 
near-infrared fluxes in this work are uncertain at the ~ 50% level. Figures 6-9 present line 
emission contour maps superposed on the continuum images at the same wavelength/grating 
setting, whereas Figures 10 and 11 display the integrated spectra. A few of the spectral 
images exhibit strong residuals due to non-uniform illumination during the spectral flat-field 
procedure. 



3.2. Ha Imaging 

Ha images of the galaxies in this sample were obtained in several observing runs. 
NGC 2388 was observed with the 1.07 m Hall telescope at Lowell Observatory, and the 
rest of the sample was observed with the Perkins 1.8 m telescope at Lowell Observatory. 
IC 10 and NGC 1569 were observed in November 1992 with a TI 800x800 CCD on loan 
to Lowell Observatory from the US Naval Observatory in conjunction with the Ohio State 
University Fabry- Perot that was used simply as a 3:1 focal reducer. The resulting scale was 
0"488 pixeU 1 . NGC 2388 was observed in December 1997 with a SITe 2048x2048 CCD, 
binned 4x4. The pixel scale was 0'. / 62 pixel" 1 . The rest of the observations with the Perkins 
telescope were made in October-December 1996 and used an 800x800 TI CCD provided 
to Lowell Observatory by the National Science Foundation. This detector was coupled to a 
4:1 focal reducer to yield a pixel scale of 0'/433 pixel -1 . NGC 6946 was observed with the 
KPNO 0.9 m telescope in May 1993, and the data were graciously provided by A. Ferguson 
(Ferguson et al. 1998). The continuum was subtracted using an image taken with a nar- 
rowband filter centered at 6409 A and width 88 A (A. Ferguson, private communication). 
January 200 KPNO 2.1 m Ha data of NGC 1569 were kindly provided by C. Martin to 
enhance our analysis (Martin, Kobulnicky, & Heckman 2002). 

All observations at Lowell Observatory used a series of redshifted Ha filters with FWHM 
of 27-32 A, and an off-band filter centered at 6440 A with a FWHM of 95 A. The off-band 
image was shifted, scaled, and subtracted from the Ha image to remove the stellar continuum 
and leave only Ha nebular emission. The images were reduced using the over-scan strip to 
remove the electronic pedestal and dome flat-field images to remove pixel-to-pixel variations. 
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Multiple images were obtained and combined to remove cosmic rays. 

The Ha emission was calibrated using the known Ha flux of NGC 604 (D. Hunter, 
unpublished spectrophotometry) and spectrophotometric standard stars with minimal Ha 
absorption features (Stone 1977; Oke & Gunn 1983). These objects agreed with each other to 
4% and the calibrations remained constant to 8% for a given filter over many observing runs. 
The Ha flux has been corrected for [N II] contamination and the position of the line within 
the bandpass, shifted for the outside temperature. We assume Ha/(Ha + [N II] )~ 0.75 
(Kennicutt 1992; Jansen et al. 2000). To account for the temperature-dependent shifting 
of the filter bandpass, we assume the bandpass shifts blueward by 0.18 A for every degree 
Celsius colder than 20° C. Fluxes for the PIFS fields-of-view are extracted from the Ha 
images. The overall absolute calibration of these fluxes is estimated at ~ 30% to account 
for PIFS-Ha aperture misalignment, unknown observing systematics, and the uncertainties 
in [N II] contamination, which exhibits a 20% dispersion in star-forming galaxies (Kennicutt 
& Kent 1983; Kennicutt 1992; Jansen et al. 2000). 



3.3. Optical Spectroscopy and Derived Abundances 

Optical long-slit spectra of six of the galaxies studied in this paper, NGC 693, NGC 1569, 
NGC 4418, NGC 7218, NGC 7771, and UGC 2855, were obtained with the Large Cassegrain 
Spectrograph (LCS) on the McDonald Observatory Harlan J. Smith reflector during observ- 
ing runs in September 1996, March 1997, and September 1997. These data are a subset of a 
larger body of optical spectra of galaxies in the ISO Key Project, which will be presented and 
discussed in more detail in a future paper (Dinerstein et al., in preparation). The 90"-long 
slit of the LCS was opened to a width of l'/5, yielding a spectral resolving power of ~800 
on the Texas Instruments 800 x 800-pixel CCD detector. Separate setups were used to cover 
the blue (~ 3680-5060 A) and red (~ 4800-7450 A) spectral regions; both set-ups included 
H/3 and the [O III] 4959, 5007 A lines. Both blue and red data were obtained for all of 
the galaxies except NGC 1569 and UGC 2855, for which only the red spectra were taken. 
Integration times are given in Table 3. Standard IRAF procedures were used to reduce the 
data, extract and calibrate apertures corresponding to selected spatial regions, and measure 
line strengths. 

The physical regions sampled by the optical spectra vary from galaxy to galaxy, and do 
not always correspond to the apertures used for the near-infrared spectroscopy. The optical 
slit was usually rotated to the position angle of the galaxy's major axis, but in some cases 
we then offset onto emission peaks in the Ha images described in § 3.2. For NGC 693, 
NGC 4418, and NGC 7771, a single spectrum representing the integrated light of the disk 
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was extracted, while for NGC 1569, NGC 7218, and UGC 2855 we extracted spectra from 
smaller apertures centered on bright disk H II regions. 

The primary quantities we derived from the optical spectra are tabulated in Table 3: 
c, the logarithmic extinction at H/5 as derived from the ratio of Ha to H/3 (which we 
used to correct the optical line ratios for reddening); the measured H/3 flux in the extrac- 
tion aperture (uncorrected for extinction); and several bright-line indices used as empiri- 
cal metallicity indicators. These include: R 23 = ([O III] 4959, 5007 + [O II] 3727) /H/?; 
R 3 = ([O III] 4959, 5007)/H/3; 32 = ([O III] 4959, 5007)/[O II] 3727; and N/O = [N II] 6584/[0 II] 3727. 
The [S II] 6717, 6731 A intensity ratios generally indicate low values for the electron density 
n e , log n e ~ 2.0 ±0.3 (cm" 3 ). 

For the three galaxies for which we obtained blue data, NGC 693, NGC 7218, and 
NGC 7771, we used the R23 calibration of Edmunds & Pagel (1984) to obtain the (O/H) 
abundances given in Table 1. If we had used the calibrations of McGaugh (1991, 1994) 
for intermediate ionization parameters, the inferred (O/H) values would be systematically 
lower by —0.1 dex, but the galaxies would retain the same order and differences in abun- 
dance values. Note the contrast between NGC 693 and NGC 7771, which have similar line 
fluxes and line-to-continuum contrast (in integrated light) but strikingly different ionic line 
ratios (Table 3), which indicate that (O/H) and (N/O) are about 0.4—0.5 dex higher in 
NGC 7771 than in NGC 693. These two galaxies give an excellent illlustration of both the 
excitation-metallicity effect which is the basis of the bright-line method, and the well-known 
correlation between metallicity and total mass or luminosity of disk galaxies (see Garnett 
2002 for a recent discussion). NGC 7218 is intermediate in abundance between NGC 693 
and NGC 7771. Since we have only red spectra for UGC 2855, we assume that it has the 
same (O/H) as NGC 693 based on the similarity of their R 3 indices. On the other hand, 
NGC 1569 has extremely strong [O III] lines. Taking R 3 as a lower limit to R 23 , we find an 
(O/H) value of about 8.2, consistent with the result of Kobulnicky & Skillman (1997). 



4. RESULTS 

4.1. Fluxes, Equivalent Widths, and Starburst Ages 

Near-infrared emission lines are securely detected in 84% (37/45) of the lines and regions 
targeted (Table 4). In general, Pa/3 is the strongest of the observed lines, and H 2 (2.122 /xm) 
the weakest (no significant H 2 emission is observed for any H II region in the sample). This 
is also true for the extinction-corrected fluxes (§ 4.3). The most notable exception to this 
rule is provided by NGC 4418, which has a surprisingly large H 2 flux, brighter in fact than 
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the other three near-infrared lines observed combined. 

The equivalent widths are generally smaller for the nuclei than for the H II regions. For 
example, the uncertainty-weighted averages for Pa/3 and Br7 are 2.9 ± 1.5 A and 5.3 ± 0.9 A 
for nuclei and 40. ± 15 A and 24 ± 6 A for H II regions. This is no surprise since the nuclear 
regions have comparatively much stronger near-infrared continuum emission from older stars. 
Ho, Fillipenko, & Sargent (1997) found a similar result in their study of H II nuclei and disk 
H II regions. 

The Pa/5 and Br7 equivalent widths are essentially ratios of line emission from H II 
regions around young stars to continuum emission from older stars. If the conditions of either 
an instantaneous starburst or recent constant star formation are met, then these equivalent 
widths are excellent age indicators. From a modeler's perspective, one needs to know the 
star formation history for at least 1 Gyr; if there are variations in the star formation rate 
with time, the continuum is difficult to constrain. On the other hand, if the recent starburst 
is dominant and has an age less than 30 Myr, its red supergiant population is often so 
strong that the older population does not significantly affect the inferred starburst age. For 
example, in the 10 6 M Q instantaneous starburst model of Starburst99 (Figure 51 of Leitherer 
et al. 1999), the peak in the K band luminosity as a function of time (for Z > 0AZ Q ) stems 
from the burst's red supergiant population; this peak at 10 7 yr characterizes the starburst. 
Compare this luminosity to the luminosity generated by a galaxy-wide continuous 1 M© yr^ 1 
star formation that occurs for, say, 10 9 yr (see their Figure 52). Considering that only a 
small portion of this underlying global luminosity is projected to spatially overlap with an 
individual H II region, Br7 equivalent widths for an H II region can indeed be relatively 
unaffected by an older population. 

Assuming that the Br7 equivalent width can be related to the effective star formation 
age, the ages derived for an instantaneous starburst (Leitherer et al. 1999) range from a 
few million years for the youngest H II regions to longer than 10 Myr for the oldest nuclear 
regions (Table 5). The H II regions tend to have more recent starbursts than nuclei. If nuclei 
in general are not defined by a single starbursting population, and instead are comprised 
of star formation episodes that span a range of ages, then it is expected that their inferred 
ages would be older than an isolated starburst. Interestingly, excluding the unusual system 
NGC 4418, there is perhaps a trend in the starburst age with mid-infrared color. Systems 
with lower f u (6.75fim)/f u (15fim), indicating more active star formation, tend to be younger 
(see Figure 12). A similar trend was previously reported for a study of galaxy circumnuclear 
regions (Roussel et al. 2001). Finally, readers familiar with post-starburst NGC 1569 may be 
surprised by its apparently young age. We emphasize that our observations centered on two 
H II regions, not the famous super-star clusters. Moreover, the differing spatial distributions 
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in the continuum and line (§ 4.2) can lead to ill-defined equivalent widths; we have extracted 
an average equivalent width, based on the entire fluxes within the PIFS field of view. 

4.2. Spatial Distribution 

For the nuclear regions where line emission is detected, the continuum and line emission 
are generally co-spatial (have similar morphology). On the other hand, there are clear spatial 
offsets between the line emission from the ionized gas and the continuum emission in the 
disk H II regions (see the maps in Figure 6-9). 

The similar distributions of the line (gas) and continuum emission (stars) in the galaxy 
nuclei may arise naturally, as a result of interstellar material settling to the center of the 
galaxy's gravitational potential well, confining the peak density of various components to a 
narrow region. Conversely, the H II regions do not provide the same kind of gravitational 
environment, so that the distribution of ionized gas is affected by a number of other factors. 
For example, interstellar material can be blown away from a region that is actively forming 
stars: the Ha map of NGC 1569 exhibits a hole centered on the two bright super starclusters, 
and the distinct "S" shape to the near- infrared line emission maps in NGC 1569 NW is also 
apparent in Ha at the same location (Martin, Kobulnicky, & Heckman 2002; the emission 
line contours for NGC 1569 NW correspond to GMC 3 of Taylor et al. 1999). 

Alternatively, there may simply be two or more populations of different ages, as in 
the Smith et al. (1999) study of the NGC 7771 circumnuclear star-forming ring and the 
Alonso-Herrero et al. (2003) study of starbursts M 82 and NGC 253. In NGC 7771, Smith 
et al. (1999) find that the 1.6 kpc diameter 10 distinct sources in both 

the radio and near-infrared continua. However, the radio and near-infrared clumps do not 
overlap, suggesting at least two different formation timescales. Most of the ring is visible in 
the PIFS field of view, and there is evidence for at least three distinct sources in our Pa/3 
map. Alonso-Herrero et al. (2003) find a continuum- line spatial mismatch for much of M 82 
and NGC 253 in Paa and [Fe II] 1.644 /zm. In addition, given that the majority of radio 
supernova remnants and compact [Fe II] 1.644 /xm sources are not spatially aligned, these 
authors suggest that there is a young and old population of supernova remnants. 

4.3. Extinction 

We derived the extinction for each region from the ratios of Ha (from the narrow-band 
imaging), Pa/5, and Br7. The values listed in Table 5 are derived by least-squares fitting 
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a screen model to the observed line ratios Ha/Br7 and Pa/3/Br7 (the two line ratios are 
fit simultaneously) and assuming the Li & Draine (2001) extinction curve. Uncertainties 
on the derived extinctions are based on the observational uncertainties and the predicted 
range in the intrinsic line ratios: Ha/Br7= 97.5 ±5.5 and Pa/3/Br7= 5.75 ±0.15 for Case B 
recombination; 100 < n c (cm~ 3 ) < 10,000; and 5,000 < T C (K) < 10,000 (Storey & Hummer 
1995). 

A small difference is found between the nuclear and H II region extinctions. The averages 
for the sample are Ay 1.4±0.7 mag for the H II regions and 2.7±0.6 mag for the nuclei. 
Our results are consistent with previous determinations. For example, Reunanen et al. 
(2000) find A y (int±MW)= 2.8 mag for the central ll"xll" of NGC 7771, and Davies, 
Alonso-Herrero, & Ward (1997) obtain Ay = 3.0 mag for the more obscured optical nucleus 
(their "region 2"). Engelbracht et al. (1996) find that the core of NGC 6946 is extincted 
by Ay(int±MW)= 4.3 mag for a foreground screen model. The internal extinction for 
NGC 1569 SE and NGC 1569 NW has been estimated by Kobulnicky & Skillman (1997) 
to be ~0.2 mag (their regions C15 and A12, respectively). The extinction for NGC 4418 is 
unconstrained by these data, since the fluxes for both Pa/3 and Br7 are upper limits. 

The range of extinctions found by Ho, Filippenko, & Sargent (1997) for a sample of over 
200 "H II nuclei" is < A v < 3. The results of Kennicutt, Keel, & Blaha (1989) for H II 
and starburst nuclei are consistent with this result. This range of extinction for optically- 
selected samples is a bit smaller than the average nuclear extinction derived for our sample 
((Ay(int)) ~3 mag), with NGC 693 being significantly higher than the average. In fact, in 
most cases the extinction values are larger than those indicated by our long-slit-derived ratio 
of H/3/Ha: (note: for a ratio of total to selective extinction of R = Ay j ' E(B — V) = 3.1, 
A v = 2.1 c). This is to be expected in regions where dust is mixed with gas, since the optical 
line emission is seen selectively from the regions of lowest extinction. As Kennicutt (1998) 
points out, infrared-selected samples reveal much more dramatic nuclear star formation, 
and V band obscuration can easily exceed several magnitudes in the cores of star-forming 
galaxies; we are able to probe deeper at near-infrared wavelengths than optical studies, and 
the more deeply embedded sources simply are more obscured. 

4.4. Star Formation and The Interstellar Radiation Field 

The number of hydrogen-ionizing photons is derived from an average of the (extinction- 
corrected) Br7 and Pa/3 luminosities 

A^H )^ 1 ] = i{7.63L Br7 [10 20 W\ + 1.40L Pa/3 [10 20 W]} (1) 
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(Leitherer & Heckman 1995). If either the Br7 or Pa/3 flux is an upper limit, then only 
the securely detected flux is used. Table 5 also expresses this number in terms of the 
equivalent number of 08.5 V stars, using a conversion from Vacca, Garmany, & Shull (1996). 
The average projected density of ionizing photons for this sample, computed from a ratio 
of the number of ionizing photons to the size of the emitting region, spans two orders of 
magnitude, from log JV(H )[s _1 pc~ 2 ] = 46.4 for NGC 7218 to 48.2 for NGC 1569 NW 
(the uncertainty on these numbers is 0.2 dex). A strong and understandable correlation is 
obtained between the mid-infrared color /^(6.75//m) / f u (15fj,m) and ionizing photon density: 
lower / I/ (6.75/xm)// J/ (15/xm), an indicator of more active star formation, corresponds to more 
intense radiation environments (Figure 12). Forster Schreiber et al. (in preparation) find 
a similar result for disk and nuclear regions in a sample of nearby galaxies. Perhaps not 
surprisingly, the one outlier in the observed trend is NGC 4418, the most unusual target in 
our sample (see § 2.6, 4.1, and 4.3). This object may falsely appear to be lacking in ionizing 
photons either because of extreme extinction, or because its H II regions are very young and 
still rich in dust which absorbs a large fraction of the ionizing photons. Finally, with the 
exception of the H II region dubbed "NGC 6946 NE," the H II regions exhibit the most 
intense radiation fields and lowest / iy (6.75/xm)// I/ (15/xm) ratios. 

The combination of N(H°) and the H 2 flux can place constraints on the dominant exci- 
tation mechanism for H 2 . Based on the results of Black & van Dishoeck (1987), I\jv/Ih 2 ~ 
10 5 — 10 8 for H 2 fluorescence triggered by ultraviolet radiation (H 2 fluorescence, however, can 
initially be much more significant compared to its final equilibrium value, e.g. Goldshmidt & 
Sternberg 1995 and Hollenbach & Natta 1995). For our sample, the observed ratio is of order 
10 4 , suggesting that the dominant excitation mechanism is shocks rather than fluorescence. 

4.5. Near-Infrared Line Ratio Diagnostics 

The [Fe II] /Pa/5 line ratio can be a useful diagnostic of the star formation history since 
[Fe II] emission is thought to arise from supernova activity. Accordingly, [Fe II] is an indicator 
of past star formation, whereas Pa/3 signifies current star formation (e.g. Moorwood & Oliva 
1988; Calzetti 1997). In a sense [Fe II] /Pa/3 is a ratio of non-thermal to thermal activity 
(however, both thermal and non-thermal mechanisms can be attributed to supernovae and 
their remnants, e.g. Sankrit et al. 1998; Rho et al. 2003). While other lines also trace 
the current star formation, Pa/5 is particularly useful since [Fe II] 1.257 /xm/Pa/3 is relatively 
insensitive to extinction. Infrared-bright galaxies that are predominantly powered by star 
formation should thus have lower [Fe II] /Pa/3 ratios than AGN-dominated systems (e.g. 
Larkin et al. 1998). 
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H 2 (2.122 / um)-to-Br7 is another ratio that in principle is unbiased by extinction. Its 
interpretation, however, is more difficult. H 2 , in contrast to [Fe II], can be triggered by both 
thermal and non-thermal processes. As described by Calzetti (1997), molecular hydrogen 
on the surface of molecular clouds can be collisionally excited in gas heated by slow shocks, 
ultraviolet radiation, and X-ray radiation, or it can result from fluorescence following ul- 
traviolet absorption. Empirically, starburst galaxies and H II regions exhibit lower H2/Br7 
ratios (Larkin et al. 1998). 

Inspection of Figure 13 confirms that these line ratios help to empirically distinguish 
between galaxies powered by different processes. Data from the literature show there is an 
approximate trend proceeding from the upper right to the lower left of increasing star for- 
mation activity and decreasing AGN activity (or, in the absence of AGN activity, decreasing 
star formation activity occurring in brief, intense and embedded bursts with shocks). In 
fact, all three H II regions for which we measured both line ratios reside in the lower left, 
and the one source in our sample suggested to harbor AGN activity, NGC 4418 (Roche et 
al. 1986; Baan, Salzer, & LeWinter 1998; Spoon et al. 2001), is the upper-rightmost data 
point in our sample. The nuclear observations for the remainder of our sample lie between 
these extremes. It should be noted, however, that systems with low metallicity might be 
expected to have low iron-to-hydrogen ratios regardless of the source that powers their in- 
frared emission; the location of IC 10 and NGC 1569 in Figure 13 must reflect in part their 
low metallicity. 



4.6. Singly-Ionized Iron-to-Hydrogen Mass Ratio 

Following Forster Schreiber et al. (2001), the gas-phase abundance of Fe 1+ can be 
estimated assuming a single-layer model of uniform density and temperature: 

L([FeII])j(H) 

n(H) L(H)j([FeII]) 1 ] 

where L is the observed line luminosity, n the number density, and j the line emissivity. 
Br7 is the reference H line in this work. The emissivities are computed in the low density 
limit 11 assuming T c = 5000 K (Storey & Hummer 1995; Hamann & Ferland 1999). The Fe 1+ 
abundances are listed in Table 5. For comparison, the total gas-phase iron abundance in the 



11 The critical densities of the observed Fe transitions are of order 10 4 cm~ 3 (Hamann & Ferland 1999; 
Quinet, Le Dourneuf, & Zeippen 1996), so we assume collisional de-excitation is unimportant when computing 
line emissivities. 
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Orion Nebula is estimated to be of order Fe/H~ 10 6 (see Rubin et al. 1997 and references 
therein) . 

4.7. NGC 4418 

One particular object exhibits extreme properties in essentially every category. NGC 4418 
is isolated in plots of near-infrared line ratios, ionizing photon density and starburst age. 
Unlike every other source, it was undetected in Pa/3 and Br7, and H 2 is the brightest of the 
four transitions sampled. It is possible that NGC 4418 deviates from the overall age-infrared 
color trend due to a more complicated star formation history. Some claim the galaxy harbors 
a heavily obscured AGN (Roche et al. 1986; Spoon et al. 2001), while others suggest it is 
a pre-starburst galaxy. Roussel et al. (2003), for example, use new multi- wavelength radio 
data to argue that galaxies like NGC 4418 may be experiencing the very onset of a starburst 
after a long time of quiescence, giving them a distinctive infrared emission signature in both 
continuum and lines. Other studies have shown NGC 4418 to likewise be peculiar (§ 2.6). 
Though our data do not have the spectral and spatial resolution to resolve this debate, they 
similarly point to an object that is not currently powered by star formation. 

5. SUMMARY AND DISCUSSION 

Near-infrared spectroscopy over two spatial dimensions has been obtained for the first 
time for nuclear and extranuclear star-forming regions in a diverse sample of nine nearby 
galaxies. The sample was previously observed with ISO as part of the Key Project on the 
Interstellar Medium of Normal Galaxies. Mid-infrared and other ancillary data are used 
in conjunction with continuum and line maps at Pa/3, Br7, [Fe II], and H 2 wavelengths to 
explore the physical processes underlying star formation. 

In this small sample, the H II regions differ from the galaxy nuclei according to nearly 
every metric. Nuclear observations typically show spatial coincidence of the continuum and 
line emission, unlike that seen for star-forming extranuclear regions. In addition, the star- 
forming episodes characterizing the H II regions appear to be younger and more intense 
than in the nuclei. One hypothesis is that the density of material in nuclear regions supplies 
continual, on-going episodes of star formation, while extranuclear star-forming regions are 
more transient and dependent upon random overdensities in the interstellar distribution. 
Thus if measures of a single star-formation age are made, such as that inferred from Br7 
equivalent widths, on-going nuclear star formation would naturally lead to longer timescales 
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than that seen for individual starbursting episodes in H II regions. Moreover, continual 
nuclear star formation means both the near-infrared continuum and line transitions would 
be strong at the same locations, even if they represent stellar populations of different ages. 
Alternatively, though there may be cycling in the nuclear star formation history, the strong 
central density gradient should result in spatially coincident continuum and line emission. 
Though the H II regions appear to differ from the nuclei by most measures, this result is 
perhaps not surprising nor profound, especially given our limited sample of targets. 

Near-infrared line ratios that distinguish between AGN and star-forming activity are 
consistent with the H II regions being dominated by star formation whereas the nuclear 
regions reflect a mixture of star formation and a secondary ionization component (in the 
absence of AGN activity, the secondary component could be due to embedded and intense 
starburst-induced shocks). Perhaps the nuclei are not exclusively powered by normal star 
formation, echoing the suggestion of Terlevich et al. (1992) that AGN activity is the natural 
consequence of the final stages of intense nuclear star formation. Ho, Filippenko, & Sargent 
(1997) find that about half (42%) of the bright spirals have nuclear emission spectra with 
H II region-like line ratios; the percentage is a steep function of global morphology, with up 
to 80% of the later types exhibiting H II nuclei. Kennicutt, Keel, & Blaha (1989) find that 
approximately "half of the H II region and starburst nuclei show evidence for a secondary 
ionization component, either an active galactic nucleus or large-scale shocks ... [supporting] 
the idea that many spiral nuclei are composite in nature, with a central LINER or Seyfert- 
like nucleus surrounded by star-forming regions." Consistent with this interpretation, our 
nuclear data do not exhibit extreme line ratios indicative of pure star formation or AGN 
dominance, but instead suggest a combination of the two processes. 

One of the secondary goals of this work was to search for any trends in the sample as 
a function of metallicity. No obvious correlations are seen with respect to the oxygen abun- 
dances. However, on a purely empirical basis the H II regions show lower gas-phase Fe 1+ 
abundances than nuclei by an order of magnitude. Total gas-phase Fe abundances cannot 
be computed with these data since only the singly-ionized species has been observed. These 
data are additionally limited since Fe 1+ abundances are difficult to interpret. Not only do 
Fe 1+ abundances depend on the overall amount of Fe produced, Fe 1+ abundances also depend 
on the ionizing strength of the interstellar radiation field and how much Fe has been depleted 
onto, or released from, dust grains. These complications are likely why the abundances do 
not correspond in any way to the derived metallicities. Nevertheless, there is an interesting 
contrast in the Fe 1+ /H abundance for the H II-like regions and the nuclei. The H II-like 
regions show low abundances in comparison to the nuclei. Lack of abundance data from ad- 
ditional ionic species leaves many questions unanswered. For example, do environments with 
harder radiation fields like H II regions exhibit relatively higher Fe 2+ and Fe 3+ abundances, 
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helping to explain their relatively low Fe 1+ abundances? Models of Orion Nebula data indi- 
rectly support this scenario, with ionization fractions of Fe 1+ /Fe 2+ /Fe 3+ =0. 0128/0. 244/0. 744 
(Baldwin et al. 1991) and 0.0529/0.414/0.533 (Rubin et al. 1991). To directly compare these 
fractions in a relatively soft and a relatively hard radiation environment, we utilize the NEB- 
ULA software (Rubin et al. 1991 as updated recently Rodriguez & Rubin 2003). Keeping 
the metallicity (Orion-like) and the number of ionizing photons (10 50 s _1 ) the same from 
simulation to simulation, we find a relatively small change in the Fe 1+ abundance but a clear 
shift to higher Fe 3+ abundances, mostly at the expense of Fe 2+ , for harder radiation fields 
(see Table 6). 

Perhaps the most interesting result is that portrayed in the bottom panel of Figure 12, 
the relation between ionizing photon density and mid-infrared color. Relatively intense 
radiation fields result in lower f u (6.75/im)/ f u (15/im), likely via two different effects: reduced 
emission from polycyclic aromatic hydrocarbons and a steepening mid-infrared continuum 
at 15 fim (e.g. Dale et al. 1999; Dale et al. 2001; Helou et al. 2001). Intense radiation 
fields could result from higher photon densities or harder radiation (Roussel et al. 2001). 
To explore this issue, we turn to 32 , the ratio of [O III]5007+4959 to [O II]3727 flux. 12 
Three de-reddened measurements are listed in Table 3, and four de-reddened values are 
drawn from elsewhere: log0 32 = 0.53, -1.04,0.88,0.51 for IC 10, the nucleus of NGC 6946, 
NGC 1569 NW, and NGC 1569 SE, respectively (Lequeux et al. 1979; Heckman, Crane, & 
Balick 1980 13 ; C. Kobulnicky, private communication). Once again, the H II regions behave 
differently than the nuclei, log0 32i Hii > 0.5 and log0 32iimc ^ —0.5, with the H II regions 
exhibiting harder radiation fields. This 32 segregation according to environment implies 
that mid-infrared colors are in part determined by the hardness of the radiation field. 

We owe gratitude to several who unselfishly contributed: J. Larkin for near-infrared line 
ratio data; A. Ferguson for an Ha image of NGC 6946; C. Kobulnicky, C. Martin and B. 
Buckalew for NGC 1569 Ha images; M. Regan for preliminary investigations of CO data; 
D.G. Hummer and H.W.W. Spoon for their program that extracts intensity ratios of hydro- 
genic recombination lines; and C. Leitherer and A. Bianchini for enlightening discussions. 
We also would like to acknowledge the Palomar staff and B.T. Soifer who helped make the 
near-infrared observations possible. Suggestions by the referee improved this work. This 
publication makes use of data products from the Two Micron All Sky Survey, which is a 



12 The parameter O32 is not the best discriminator for radiation hardness for systems with 12 + log(0/H) > 
8.7, like NGC 7771 and NGC 6946 (see Figure 8 of Kobulnicky, Kcnnicutt, & Pizagno 1999). 

13 The value for NGC 6946 was derived assuming that the [O III]4959 flux is 1/3 of the quoted [O III]5007 
flux. 
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-0.48 


8.1 


8.17 


NGC 0693 


10: sp 


015031.0 


+060842 


22.1 


1567 


+0.13 


9.8 


8.4 


UGC 2855 


SB(s)cd II-III 


034822.6 


+700757 


18.7 


1202 


+0.31 


10.4 


8.4 


NGC 1569 


IBm; Sbrst 


043049.0 


+645053 


2.5 


-104 


-0.66 


8.7 


8.19 


NGC 2388 


SA(s)b: pec 


072853.5 


+334905 


54.8 


4134 


+ 1.11 


10.9 




NGC 4418 


(R')SAB(s)a 


122654.6 


-005240 


27.3 


2179 


+1.12 


10.6 




NGC 6946 


SAB(rs)cd 


203452.3 


+600914 


4.5 


48 


-0.34 


9.8 


9.28 b 


NGC 7218 


SB(r)c 


221011.7 


-163936 


23.8 


1662 


-0.17 


9.7 


8.6 


NGC 7771 


SB(s)a 


235124.8 


+200642 


60.0 


4277 


+0.41 


11.1 


9.1 



a Dalc ct al. 2000. 

b Assuming the metallicity trend with radius found by Pilyugin & Ferrini 1998, the targeted H II region has 
logO/H + 12 = 9.08. 

c Metallicity references: IC 10: Lequeux et al. 1979; NGC 1569: Kobulnicky & Skillman 1997; NGC 6946: Pilyugin 
& Ferrini 1998; NGC 693, UGC 2855, NGC 7218, NGC 7771: this work 
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Table 2. Palomar 200 inch Telescope Observations 



Date 


Oalaxv a 


R A Decl 


Lines* 3 


\\ band 


Skv 






fJ2000) 




SeeinEf 


Conditions 


2000 Sep 16/17 


NGC 6946 


203452.3 +600914 


[FeIl],Pa/5\H 2 ,Br 7 


0"9 


Photometric 




NGC 1569 NW 


043047.3 +645101 


Pa/5 


0"9 


Photometric 


2000 Sep 17/18 


NGC 6946 NE 


203506.0 +601057 


[FeIl],Pa/5\H 2 ,Br 7 


0"6 


Thin cirrus 




NGC 1569 NW 


043047.3 +645101 


[Fe II] 


0"6 


Photometric 


2000 Sep 18/19 


NGC 6946 N 


203452.3 +600918 


[FeIl],Pa/5\H 2 ,Br 7 


0"7 


Thin cirrus 


2000 Oct 13/14 


NGC 6946 S 


203452.3 +600910 


[Fell], Pa/3 


0T8 


Thin cirrus 




NGC 1569 NW 


043047.3 +645101 


[Fell], H 2 ,Br 7 


0"8 


Thin cirrus 


2000 Oct 14/15 


NGC 6946 S 


203452.3 +600910 


[Fell], H 2 ,Br 7 


0"7 


Photometric 




NGC 6946 


203452.3 +600914 


[Fe II] 


0"7 


Photometric 




NGC 1569 SE 


043051.7 +645048 


Br 7 


0"6 


Photometric 


2001 Jan 06/07 


NGC 1569 SE 


043051.7 +645048 


Pa/5 


0"5 


Thin cirrus 




NGC 2388 


072853.5 +334908 


[FeIl],Pa/5\H 2 ,Br 7 


0"5 


Thin cirrus 




NGC 4418 


122654.6 -005240 


[FeIl],Pa/5\H 2 ,Br 7 


0"6 


Thin cirrus 


2001 Oct 03/04 


NGC 7218 


221011.7 -163939 


[FeIl],Pa/5\H 2 ,Br 7 


0"8 


Photometric 




UGC 2855 


034820.4 +700758 


[Fell], Pa/5, Br 7 


0"9 


Photometric 


2001 Oct 04/05 


NGC 7771 


235124.8 +200642 


[FeIl],Pa/5\H 2 ,Br 7 


0"9 


Photometric 




IC 10 


002027.2 +591743 


Pa/5,H 2 ,Br 7 


1"2 


Photometric 


2001 Oct 05/06 


IC 10 


002027.2 +591743 


[Fe II] 


1"6 C 


Thin cirrus 




NGC 0693 


015030.9 +060842 


[FeII],Pa/5,H 2 ,Br 7 


0"7 


Thin cirrus 


2001 Nov 05/06 


UGC 2855 


034820.4 +700758 


H 2 


0"8 


Photometric 



a Bold lettering indicates the initial observation(s) for a source. If a source was targeted more than once, additional 
entries are not indicated in bold. 

b The [Fe II] line observed was either [Fe II] 1.257 or [Fe II] 1.644 fim, or both. 
C H band measurement. 
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Table 3. Parameters Derived from Long-slit Optical Spectra 



Galaxy 


Date 


T cxp (b,r) log(H/3) a c(Ha,H/3) b log R 23 c log R 3 d 


log 32 c 


loj 




Obs. 


(minutes) 







NGC 0693 


Sep 1996 


40 


90 


-13.6 


1.05±.10 


0.84±.02 


+0.25±.02 


-0.47±.03 


UGC 2855 


Sep 1997 




50 


-14.9 


2.0 ±.2 




-0.26±.10 




NGC 1569 


Sep 1997 




20 


-12.9 


1.1 ±.1 




+0.94±.01 




NGC 4418 


Mar 1996 


40 


45 












NGC 7218 


Sep 1996 


30 


90 


-14.3 


0.5 ±.1 


0.69±.03 


-0.02±.01 


-0.62±.04 


NGC 7771 


Sep 1996 


60 


60 


-13.8 


1.1 ±.1 


0.28±.10 


-0.51±.02 


-0.71±.12 



a Observed flux in erg cm 2 s 1 , for the selected aperture. 

b c = log 10 [/(H/3)t r uo//(H/3)obs]- For a ratio of total-to-selective extinction of R — A V /E(B - V) = 3.1, A v = 2.1c. 

c The commonly-used bright-line index, R 23 = QO III] 5007+4959 + [O II] 3727)/H/3; no value is listed for objects 
without blue data. 

d R 3 = QO III] 5007+4959)/H/3. 

c 32 = ([O III] 5007+4959)/ [Oil] 3727. 

f N/0 = [Nil] 6584/[OII] 3727. 
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Table 4. Line Fluxes 3, and Equivalent Widths b 



^aiaxy 


i ai get 


ft™ 

net 

II (~\^\(~\ It Y"fl 

U.UOU /JillL 


rrp tt| 
ire hi 

1 9^7 nm 


Pa R 
1 981 8 nm 


\T?a TTl 

[re iij 

l.U^^ yU.111 


XJ 

112 


nvy 

J a II Wi 

Z.1UOO /Jilll 


IC 10 


HII 


1400 




200 


<1.5 


<0.7 


17 










190± 7 


<5.0 


<2.5 


59± 2 


NGC 0693 


nucleus 


96 




3.9 


4.5 


3.7 


5.4 










1.1±0.1 


1.8±0.1 


2.7±0.2 


3.7±0.2 


UGC 2855 


nucleus 


34 




90. 


41 


8.9 


31 










32± 2 


19± 1 


3.6±0.2 


24± 2 


NGC 1569 NW 


HII 


1500 


6.6 


420 


9.1 


6.4 


100 








1.8±0.1 


100± 2 


3.6±0.2 


5.2±0.3 


82± 3 


NGC 1569 SE 


HII 


400 




55 






14 










30.8±0.7 






21.9±0.4 


NGC 2388 


nucleus 


180 


18 


68 


26 


6.7 


12 








3.5±0.1 


18.5±0.7 


13.9±0.8 


2.8±0.6 


5.3±0.8 


NGC 4418 


nucleus 


5.5 




<0.5 


0.7 


16 


<1.1 










<0.2 


0.5±0.1 


13.9±0.5 


<0.9 


NGC 6946 


nucleus 


290 


88 


180 


75 


61 


64 








2.8±0.1 


9.0±0.2 


10.5±0.3 


6.5±0.2 


7.0±0.2 


NGC 6946 NE 


HII 


50 


<1.6 


22 




<2.1 


8.5 








<13 


41± 4 




<23 


140±40 


NGC 7218 


nucleus 


66 




8.0 


1.9 


<0.5 


1.2 










4.8±0.4 


1.8±0.2 


<0.7 


2.4±0.7 


NGC 7771 


nucleus 


160 




40 


12 


26 


15 










7.3±0.3 


9.7±0.8 


16± 2 


5.8±0.4 


A x /A v c 




0.7834 


0.2903 


0.2815 


0.1907 


0.1247 


0.1204 



a Fluxes corrected for atmospheric extinction are given in the first row for each galaxy and are in units of 
10~ 18 W toT 2 . The Ha fluxes (~ 30% uncertainty) are extracted from image cut-outs matched to the PIFS 
fields-of-view. The absolute calibration of the near-infrared fluxes is <~ 50% (see text). 

b Equivalent widths are listed in the second row for each galaxy are in units of A. 



c Dcrivcd from Li & Drainc 2001. 
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Table 5. Results 



Galaxy 


larget 


i /„(6.75«m) a 
1U & /,(15 M m) 


A y (MW,int) b log N(H°) 


1 ^ ATC 


a d 

Age SB 


TT 1 + /TTi 

be^/ri 








I 5 J 




(10 7 yr) 


(10 7 ) 




H TT 

11 11 




2.39, 0.0±2.7 


50.5 


1 77 


n r 




NGC 0693 


nucleus 


-0.07 


0.17,11. 7±4.4 


53.1 


4.35 


0.8 


25 


UGC 2855 


nucleus 


-0.19 


2.45, 3.7±0.8 


53.4 


4.72 


0.7 


26 


NGC 1569 NW 


HII 


-0.93 


2.32, 0.3±0.9 


51.9 


3.23 


0.4 


0.82 


NGC 1569 SE 


HII 


-0.69 


2.32, 0.0±3.8 


51.1 


2.41 


0.5 




NGC 2388 


nucleus 


-0.34 


0.19, 2.3±0.8 


53.8 


5.05 


0.7 


19 


NGC 4418 


nucleus 


-0.79 


0.08, ••• f 


<52.7 


<3.98 


>0.9 


>23 


NGC 6946 


nucleus 


-0.37 


1.13, 3.1±0.9 


52.3 


3.62 


0.7 


14 


NGC 6946 NE 


H II 


+0.02 


1.13, 2.8±0.9 


51.4 


2.69 


0.5 


<1.1 


NGC 7218 


nucleus 


-0.07 


0.11, 0.6±0.9 


51.9 


3.19 


0.8 


20. 


NGC 7771 


nucleus 


-0.17 


0.24, 3.0±0.9 


53.8 


5.12 


0.7 


12 



a This mid-infrared ratio stems from flux densities extracted over a 9" diameter circular aperture, using images 
from Dale et al. 2000 that are smoothed to ~9" resolution. 

b The foreground extinction due to the Milky Way is listed first, followed by the internal extinction. The foreground 
extinctions are compiled from the Schlegel, Finkbcincr & Davis 1998 values listed in NED, except for low Galactic 
latitude IC 10, for which the value is uncertain. For IC 10, we use the value derived by Massey & Armandroff 1995 
and suggested by Hunter 2001. 

c Derived from \ogN(H°) (e.g. Vacca, Garmany, & Shull 1996). 

d Compiled from the instantaneous starburst curves in Leitherer & Hcckman 1995 for the Br7 equivalent widths. 

e The Fe 1+ gas-phase abundance is computed using Equation 2 and assuming T e = 5000 K for the following line 
emissivitics in 10~ 22 W m~ 3 : j([Fe II]1.257 /xm)= 26.8 and j([Fe II]1.644 fj,m)= 5.94 (Hummer et al. 1993; Pradhan 
& Zhang 1993); j( Br 7) = 7.23 • 10~ 7 applicable for the range 100 < n c (c m - 3 ) < 10,000 (Storey & Hummer 1995). 
The uncertainties are estimated to be a factor of two. Three targets were observed at both 1.257 and 1.644 /xm, 
and the separately-derived Fe 1+ abundances (in 10~ 7 ) are consistent with the average at this level: 0.3 vs 1.3 for 
NGC 1569 NW; 6.5 vs 32 for NGC 2388; 8.1 vs 20 for NGC 6946. 

r 

The extinction for NGC 4418 is unconstrained; the extinction based on upper limits for Pa/3 and Br7 is 
Ay(int) = 14.4 mag. In fact, many of the results for this peculiar system should be interpreted with caution, in- 
cluding the starburst "age" (see § 4.1). 
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Fig. 1. — 2MASS K band image (5f3x5'3; J2000 coordinates) of IC 10 with two overlays: 
the PIFS field of view outlined in white and the B Band 25 mag arcsec -2 isophote in black. 



Fig. 2.— Similar to Figure 1 for NGC 693 and UGC 2855 (l.'5xi:5 2 MASS K band images). 



Fig. 3.- Similar to Figure 1 for NGC 1569 and NGC 2388 (If5xlf5 2MASS K band 
images) . 



Fig. 4.— Similar to Figure 1 for NGC 4418 and NGC 6946 and 6f2x6^2 2MASS 

K band images). 



Fig. 5.- Similar to Figure 1 for NGC 7218 and NGC 7771 (V.hxl'.h 2MASS K band 
images) . 



Table 6. Simulated Fe Ionization Fractions Fe 1+ /Fe 2+ /Fe ; 



T cfr (K) 


n e = 10 cm 3 


n e = 100 cm 3 


35,000 


0.0793/0.755/0.166 


0.0360/0.752/0.212 


50,000 


0.0746/0.218/0.707 


0.0380/0.127/0.835 



a Thcse numbers correspond to an Orion-like metallicity and 10 50 s 
ionizing photons. 
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Fig. 6. — PIFS continuum and line maps for IC 10 (top four panels), NGC 693 (middle), 
and UGC 2855 (bottom). In each panel the continuum image is shown in greyscale and the 
line emission is portrayed as contours. The contours are plotted at levels of 2, 3, 5, 7, 10, 
14, 20, and 28a. 
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Fig. 7.— Similar to Figure 6 for NGC 1569 NW (top four panels), NGC 1569 SE (middle), 
and NGC 2388 (bottom). 
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NGC 4418 
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arcseconds arcsecords 




NGC 6946 
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arcseconds arcseconds 



Fig. 8- 



Similar to Figure 6 for NGC 4418 (top four panels) and NGC 6946 (bottom). 
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NGC 6946 NE 
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arcseconds arcsecords 




NGC 7218 
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arcseconds arcseconds 
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Fig. 9.— Similar to Figure 6 for NGC 6946 NE (top four panels), NGC 7218 (middle), and 
NGC 7771 (bottom). 
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Rest Wavelength (/xm) 

1.26 1.28 1.3 1.62 1.64 1.66 1.682.1 2.12 2.14 2.16 2.18 2.2 




I I I I I I I I I I L_J I I I I I I I I I I I I 

1.26 1.28 1.3 1.64 1.66 1.68 2.12 2.14 2.16 2.18 2.2 

Observed Wavelength (/xm) 



Fig. 10.— The spectral data are plotted as the flux density surface brightness versus 
both rest and observed wavelength (top and bottom axes). The lefthand panels display 
[Fe II] 1.257 /iin and Pa/3 data, the middle panels show [Fe II] 1.644 /xm data, and the righthand 
panels give the H2 and Bvy data. Small vertical lines indicate the line wavelengths. Small 
insets show the full profiles for IC 10 Pa/3 and NGC 1569 NW Pa/5. 
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Rest Wavelength (/i,m) 

1.86 1.88 1.3 1.38 1.64 1.66 " 1.68 8.13 8.14 8.16 3.18 3.8 8.33 




1.26 1.28 1.3 1.64 1.66 1.68 2.12 2.14 2.16 2.18 2.2 



Observed Wavelength (/-im) 



Fig. 11. — Similar to Figure 10 for the remaining spectra. A small inset shows the full 
profile for NGC 6946 Pa/3. 
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Fig. 12.— Top: The starburst ages, assuming an instantaneous starburst model, as a 
function of mid-infrared. Bottom: The density of ionizing photons strongly correlates 
with mid-infrared color /y(6.75//m)// i ,(15//m), consistent with the interpretation that lower 
/ I/ (6.75/xm)// I/ (15/xm) signifies more active star formation. 
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Fig. 13. — Extinction-corrected near-infrared line ratios for this sample and for other objects 
drawn from the literature (Larkin et al. 1998; Calzetti 1997). The data for NGC 4418 have 
not been extinction corrected due to the extinction uncertainty. A sample of 11 Galactic 
compact star-forming regions spans the range 0.4 <H 2 /Br7< 2.2 (Smutko & Larkin 1999). 
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